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ABSTRACT 

We present Integral Field Unit (IFU) spectroscopy of four brightest cluster galaxies 
(BCGs) at z ~ 0.1. Three of the BCGs have close companions within a projected radius 
of 20 kpc and one has no companion within that radius. We calculate the dynamical 
masses of the BCGs and their companions to be 1.4 x 10 11 < Md yn (M Q ) < 1.5 x 10 12 . 
We estimate the probability that the companions of the BCGs are bound using the 
observed masses and velocity offsets. We show that the lowest mass companion (1:4) 
is not bound while the two nearly equal mass (1:1.45 and 1:1.25) companions are 
likely to merge with their host BCGs in 0.35 Gyr in major, dry mergers. We conclude 
that some BCGs continue to grow from major merging even at z ~ 0. We analyse the 
stellar kinematics of these systems using the \r parameter developed by the SAURON 
team. This offers a new and unique means to measure the stellar angular momentum 
of BCGs and make a direct comparison to other early-type galaxies. The BCGs and 
their companions have similar ellipticities to those of other early-type galaxies but 
are more massive. We find that not all these massive galaxies have low as one 
might expect. One of the four BCGs and the two massive companions are found to be 
fast-rotating galaxies with high angular momentum, thereby providing a new test for 
models of galaxy evolution and the formation of Intra-Cluster Light. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: 
kinematics and dynamics — galaxies: clusters: general 



1 INTRODUCTION 

Brightest Cluster Galaxies (BCGs) include the most mas- 
sive galaxies in the Universe. Models of hierarchical struc- 
ture formation naturally feature the ongoing growth of the 
most massive galaxies by mergers ( |Peebles fc Yu|1970[ ), and 
BCGs are predicted to have undergone more mergers than 
less massive galaxies (e.g. |De Lucia fc BL aizot 200f| Thus de- 
termining the merging history of BCGs is a particularly sen- 
sitive test of current formation models. As the semi-analytic 
models of |Bower et al.| ( |2006[ ); |Croton et aL] ( |2006| > show, 
massive galaxies are over-produced in N-body dark matter 
cosmological simulations, and feedback mechanisms are re- 
quired to bring the luminosity function into agreement with 
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observations. Tracing the recent assembly history of BCGs 
is vital to future development of galaxy formation models. 

Observationally the evidence for BCG growth is con- 
tradictory: Studies of the luminosities and stellar masses 
of BCGs show little evolution in mass since z ~ 1 — 1.5 
(e.g. |Brough et al.|2002| |Collins et al.|2009[ ) and their steep 
metallicity gradients are consistent with passive evolution 
since z ~ 2 (Brough et al. 2007). However, the large radii 
and low surface brightnesses of BCGs compared to normal 
elliptical galaxies are consistent with products of major, dis- 
sipationless mergers (e.g. |Oegerle fc Hoessel|[l991| |Brough| 
et al.||2005| |von der Linden et al.||2007| |Lauer et al.||200"7] r 
Their sizes and velocity dispersions may have also evolved 
faster than less-massive early- type galaxies since z ~ 0.3 
( Bernardi|2009"l although, c.f. |Stott et al.|2011[ ). While BCGs 
are frequently observed to have multiple nuclei and close 
companions (e.g. Schneider et al.|1983{ , there are only a few 
examples where the companions have spectroscopically been 
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Table 1. Properties of clusters from von der Linden et al. (2007) 
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confirmed to be bound (e.g. Tran et al. 2008 Rasmussen 
eTaLpHO] ). BCG merger histories are also strongly linked 
to the formation of Intra-Cluster Light (ICL) as it is likely 
that a fraction of any galaxy merging with the BCG ends up 
in the ICL (e.g. |Conroy et al.||2007| [Puchwein et al.|2lH0| ). 

Integral Field Unit (IFU) spectroscopy has opened a 
new parameter space in which to analyse the stellar kine- 
matics of early-type galaxies. The SAURON team have de- 
veloped a new parameter, Xr, which utilises the increased 
spatial information from IFU spectroscopy to quantify the 



observed st ellar angular momentum in galaxies (Emsellem 
et al.||2007| l: High-mass (M dyn > 2 x 10 11 Mq) early-type 
galaxies are dominated by slow rotators with low angular 
momentum, while lower mass early-type galaxies tend to 
be fast rotators with high angular momentum. While this 
parameter provides a unique way to compare the stellar 
kinematics of BCGs to those of other early-type galaxies, 
the SAURON sample has only three galaxies with Af dyn > 
6 x 10 11 Mq and only includes one BCG (M87). Observa- 
tions of BCG stellar kinematics to date have been limited to 



long-slit spectroscopy, (e.g. Kelson et al. 2002 ). The majority 
of BCGs do not show significant rotation in these observa- 
tions, but a small number of massive BCGs are observed to 



undergo some rotation about their major-axes (e.g. Brough 
et al.||2007| |Loubser et aT1|2008[ ) thus it remains unclear as 
to what angular momentum BCGs have. 

In this letter we present the first targetted observations 
of BCG stellar kinematics from IFU spectroscopy of three 
BCGs with close companions and one control BCG with no 
close companion. We use the kinematic information to un- 
dertake a detailed analysis of whether the systems are likely 
to merge. We then map the stellar kinematics of the BCGs 
and their companions and classify their angular momentum 
compared to the early-type galaxies of the SAURON sample. 

Throughout this letter we assume a Hubble constant of 
H = 70 km s _1 Mpc" 1 and Q M = 0.3, tt A = 0.7. 



2 OBSERVATIONS 

In 



von der Linden et al. (20071, a sample of 625 BCGs 



(z ^ 0.1) were selected from the C4 cluster catalogue (Miller 
et al.||2005| of the Third Data Rele ase of the Sloan Digital 
Sky Survey (SPSS; [York et al.|2000| >. We have visually iden- 
tified those BCGs with companions within ~ 10" (18 kpc 
at z ~ 0.1). Around 20 per cent of the sample have visually 
identified massive companions within ~ 10". We observed 
three BCGs with companions, plus a fourth with no com- 
panion within the same radius as a 'control' object. The 
properties of the four host clusters are described in Table [l] 
The BCGs were observed in 2008 April and May with 



VIMOS. VIMOS was used in the IFU mode with the high- 
resolution, blue grism and a spatial sampling of 0.67"/pixel. 
This gives a field-of-view of 27" x 27" and a wavelength 
range of 4200 — 6200A, with spectral resolution of 2.1 A at 
5100A (corresponding to Oi nst = 52 kms -1 ). Observations 
were made during dark time, with an average seeing of 0.7" 
(FWHM). Each BCG field was observed for 3 x 1150s. 

Initial data reduction was achieved using the VIPGI 
pipeline ( Scodeggio et al.|2005 1 . The wavelength calibration 
is accurate to 0.09A (from the 5577A sky line). Each ob- 
servation consists of three dithered exposures split into four 
quadrants. For each quadrant, a sky background spectrum 
was calculated by taking the median over pixels without 
galaxy light. The three exposures were combined with a 5a- 
clipped mean. We used the 2-dimensional adaptive spatial 



binning code of Cappellari & Copin ( 2003 1 to re-bin pixels 



to a minimum signal-to-noise ratio of 10 per pixel. 

The stellar kinematics (velocity, V, and line-of-sight ve- 
locity dispersion, a) of each galaxy were computed from 
the spectra of each bin, using a penalized pixel fitting 
scheme, (pPXF; Cappellari & Emsellem 2004 1 and the 



MILES (Medium-resolution Isaac Newton Telescope Library 
of Empirical Spectra; |Sanchez -Blazqucz et al. 2006| evolu- 
tionary stellar population templates. These templates cover 
a similar wavelength range to VIMOS and have a similar 
spectral sampling (FWHM = 1.8 A). We determined the op- 
timal penalty value of 0.2 for pPXF by requiring the max- 
imum bias in the GaussHermite parameters h3 and h4 to 
be equal to rms/3, where the rms is the scatter of the mea- 
surements obtained from Monte Carlo simulations with the 
adopted S/N and a well-resolved stellar dispersion a > 180 
km s ' 1 . The choice of penalty value does not significantly al- 
ter the result. As an external check we compare the velocity 
measured in a central aperture of 3" for each galaxy with 
those from SDSS. This gives a mean difference for the three 
BCGs and three companion galaxies with SDSS redshifts 
(TableB of cz IFU - czsdss = 37 ± 58 kms -1 . 



3 RESULTS 

3.1 Are the companions bound? 

We first examine whether the companions are bound to the 
BCGs. We estimate the probability of the systems being 
bound by considering the solid angles over which the sys- 
tems would be bound, gi ven the observed veloc ities and 

Table ^ 



masses (M dyn = 5R a a e /G; Cappellari et al.[2006 



and the Newtonian binding criterion that a two-body sys- 
tem is bound if the potential energy of the system is equal 
to or greater than the kinetic energy (e.g. |Beers et al.|1982| 
Brough et al.|2006 l 



We first analyse the 2086 system which has a relatively 
low-mass companion. The 2086 system is separated by a pro- 
jected distance of 24 kpc (2.2/Z e! BCG), 565 kms -1 and a mass 
ratio of 4:1. This system is only likely to be bound at the 
0.1 per cent level, suggesting that the low-mass companion 
is not bound to the BCG and is unlikely to merge. 

Examining the systems with massive companions: The 
1027 companion is at a projected distance of 18.1 kpc 
(2.3i? e ,BCG) and 261 kms -1 from the BCG and has a mass 
ratio of 1.45:1. This system is likely to be bound at the 61 
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Table 2. Properties of BCGs from von der Linden et al. (2007), SDSS and IFU data. Log(M* vdL ) is the stellar mass from von der Linden 
et al. (2007), measured from Petrosian magnitudes. zsDSS; SDSS an d e SDSS are the redshift, effective radius (radius containing 50 per 
cent of the Petrosian flux) and ellipticity (measured from the flux-weighted second moments, Q and U) from the SDSS. ziptj, \R e iftj, 
(T e jFTj and Log(Af(jyn,lFu) are the redshift, A, velocity dispersion and dynamical mass measurements within i? e ,SDSS from the IFU data. 
BCG 1027 has no redshift measured in SDSS. 
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Figure 1. Velocity maps scaled to 15" X 15". Left to right: Control BCG 1050, and BCGs 1027, 1066 and 2086. The centres of the BCGs 
are indicated by squares and the companions by crosses. The orientation of the images is given by the arrows. The BCGs in clusters 
1027, 1050 and 2086 do not appear to be rotating, whilst the massive BCG companions and the BCG in cluster 1066 appear to rotate. 



per cent level, with a 99.993 per cent probability that the 
companion is on a bound inbound orbit. The 1066 system is 
separated by a projected distance of 17.5 kpc (1.6i? e ,BCG) 
and 32 kins" 1 and a mass ratio of 1.25:1. It is likely to be 
bound at the 98 per cent level with a similarly high 99.989 
per cent probability that the companion is falling into the 
BCG. These probabilities suggest that these systems are 
likely to merge. The 1027 and 1066 systems are observed 
to have luminosity residuals characteristic of interactions in 
the surface brightness analysis of C4 BCGs by Liu et al. 
(2009; using the data release 2 version of the C4 cluster 
catalogue these clusters are identified as 1026 and 1055). In 
contrast, the 2086 system does not make their merger candi- 
date sample. Our detailed, kinematic analysis provides sup- 
port for their photometric method of selecting merger can- 
didates. The dynamical friction formula of Boylan-Kolchin 
et al. (2008) suggests that, depending on the eccentricity of 



the orbit, these systems would take between 0.2 and 0.35 
Gyr to merge, consistent with estimates made using a vari- 



ety of methods by Liu et al. (20091. There are no emission 



lines in our spectra so these are likely to be dry mergers. The 
companion galaxies also have very high dynamical masses, 
meaning that if they merge they will undergo nearly equal- 
mass major mergers. 



This is concrete evidence that some BCGs continue to 
grow by major, dry, mergers, even at z ~ 0. 



3.2 Stellar Kinematics 

We present two-dimensional velocity maps in Figure [l] The 
lack of a well-ordered velocity gradient across the face of 
BCGs 1027, 1050 and 2086 suggest that they do not rotate, 
as would be expected for massive early-type galaxies. In con- 
trast, BCG 1066 and the two massive companion galaxies 
appear to be rotating. The low signal-to-noise of the BCG 
2086 companion precludes detailed kinematical analysis. 

We quantify the angular momentum of these systems 
and how that relates to other early-type galaxies using the 
SAURON X R parameter ( jEmsellem et al.|2007[ ). Xr was de- 
veloped as a proxy for the observed projected stellar angular 
momentum per unit mass: Xr = (R\V\) / (RyJ\V 2 + cr 2 )). 
R is the observed distance to the centre of the galaxy, V 
the velocity, and a the line-of-sight velocity dispersion. The 
brackets correspond to a flux- weighted sky average. Xr tends 
to unity when the mean stellar rotation dominates. We esti- 
mate the error in Xr through a 10,000 iteration Monte Carlo 
simulation of the measured errors in V and a on Xr and find 
the error to be of the order of a few per cent. 

We show Xr profiles for our galaxies compared to the 
48 SAURON early-type galaxies in Figure [2] The profiles 
of BCGs 1027, 1050 and 2086 are consistent with the slow 
rotators of the SAURON sample while BCG 1066 and the 
massive companion galaxies are consistent with the fast ro- 
tators. We note that the data for BCG 1066 only reach a 
small fraction of its effective radius and this profile could 
turn over at larger radii. However, the steep profile is un- 
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Figure 2. Radial A^j profile for the 4 BCGs (thick red dashed 
lines) and 2 massive companion galaxies (thick blue lines) . The 48 
early-type galaxies of the SAURON sample are shown in grey thin 
lines with fast (slow) rotators shown by the solid (dashed) lines. 
BCGs 1027, 1050 and 2086 are consistent with the slow rotators 
of the SAURON sample whilst BCG 1066 and the companion 
galaxies arc consistent with the fast rotators. 



which, being a volume-limited sample, has limited numbers 
of high-mass galaxies ( Cappellari et al.|2010 1 . We show that 
not all massive BCGs have low angular momentum as one 
might expect. In fact, one out of our four BCGs is a fast 
rotator. We speculate here that the angular momentum of 
the BCG with the highest A_r c is related to its interaction 
with a very close, massive merging companion. However, we 
do not observe a similar increase in the 1027 system which 
also has a close, massive companion so the rotation of BCG 
1066 could be intrinsic. The two massive companion galaxies 
also show evidence for rotation. 

This investigation serves as a pilot study. In the future 
we intend to investigate a larger sample of BCGs with close 
companions to further examine whether they are bound, and 
to test the hypothesis that higher A_r c in massive galax- 
ies could be an indication of ongoing/recent merging. We 
also intend to compare directly to the high resolution (< 1 
kpc) hydrodynamical cosmological simulations of a Virgo- 
like galaxy cluster of Teyssier et al. (20101 to further under- 



stand these results from a theoretical perspective. 



likely to be due to a kinematically distinct core (KDC) as 
these are generally ~kpc sized ( jEmsellem et al.]|2007[ ) and 
would only be visible within 0.1 R e for BCG 1066. 

We take a measure of Xr at the effective radius, Xr s , or 
the largest radius in our data (c.f. [Emscllcm et al.|2007| Ta- 
ble^. BCGs 1027, 1050 and 2086 have similar X Re values. In 
contrast, BCG 1066 and the massive companions have signif- 
icantly higher values. Figure [3] illustrates where our systems 
lie with respect to the SAURON sample in terms of their el- 
lipticity (from SDSS) and dynamical mass. The BCGs have 
low ellipticities compared to the SAURON galaxies. We find 
them to have high dynamical masses showing that we are 
probing the massive end of the galaxy distribution. The low 
X Re of control BCG 1050 and BCGs 1027 and 2086 is consis- 
tent with the slow rotating SAURON galaxies while the high 
X Ra of BCG 1066 and the companion galaxies is consistent 
with the fast rotating SAURON sample. 



4 DISCUSSION 

We observe here two BCGs at z ~ 0.1 that are undergoing 
major mergers. This is evidence that some BCGs grow via 
major dissipationless merging even at z ~ 0. We also show 
that one BCG with a close companion in projection is not 
bound, demonstrating the need for kinematic analysis. 

These observations are in contrast to studies that claim 
BCGs have not increased in mass since z ~ 1 (e.g. |Whiley| 



et al.||2008| |Collins et aT|2009| |Stott et al.||2010[ ). However, 



high-resolution hydrodynamical simulations predict that in 
a major merger with a BCG, not all the companion galaxy 
merges with the BCG: 50-80 per cent of its mass will join the 
ICL (|Conroy et al.|2007| |p"uchwein et al.|2010[ ). We therefore 



infer that most of the stellar mass from late major merging 
goes into the ICL. This is consistent with observational mea- 



surements of the ICL by Gonzalez et al. ( 2005 1 



With IFU observations we open a new parameter space 
with which to test BCG formation models. Our very massive 
BCGs complement the SAURON results (in whic h there are 
only three massive galaxies Mdyn > 6 x 10 11 Mq; Emsellem 
et al.||2007| Figure jt} and the upcoming ATLAS3D results 



5 CONCLUSIONS 

We have analysed the first observations of stellar kinematics 
from Integral Field Unit (IFU) spectroscopy of four brightest 
cluster galaxies (BCGs) at z ~ 0.1. Three of the BCGs have 
companions within a projected radius of 20 kpc and one has 
no companion within that radius. 

Undertaking a detailed analysis of whether the compan- 
ions are bound to the BCGs we find new evidence that some 
BCGs grow via major dissipationless merging at z ~ 0, with 
two of the three multiple systems in our sample having a high 
probability that they are dynamically bound and likely to 
merge within ~ 0.3 Gyr. For BCGs to undergo mergers and 
also be consistent with studies that claim no evolution in the 

1 (e.g. 



Whiley et al. 2008 



stellar mass of BCGs since z 

|Stott et al.|2010| suggests that most of the stellar mass from 
the late merging goes into the Intra- Cluster Light. 

The A_r c parameter developed by the SAURON team 
offers a new and unique means to measure the projected 
stellar angular momentum per unit mass of galaxies. It en- 
ables a direct comparison of the stellar kinematics of BCGs 
to those of other early-type galaxies. The BCGs and their 
companions have similar ellipticities to those of other early- 
type galaxies in the SAURON sample but they are more 
massive. We find that not all BCGs have low angular mo- 
mentum as might be expected from an extrapolation of the 
SAURON sample. One of our three multiple-nuclei BCGs 
is a fast rotator, as are two of their similarly massive com- 
panions. These observations provide a new test for modeling 
how the most massive galaxies in the universe form. 
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Figure 3. Af} e versus (a) cllipticity, e, and (b) dynamical mass, M<j yn , for the 4 BCGs and 2 massive companion galaxies (solid points) 
with the typical uncertainty given by the error bar. The 48 early-type galaxies of the SAURON sample are divided into fast rotators 
(crosses) and slow rotators (squares). BCGs 1027, 1050 and 2086 arc consistent with the SAURON slow rotators whilst BCG 1066 and 
the companion galaxies are consistent with the fast rotators. The BCGs and their companions are rounder and more massive than the 
majority of the SAURON sample. 
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